import numpy


CBAR = 1.132000e+01
XIXX = 9.496000e+03
XIYY = 5.581400e+04
XIZZ = 6.310000e+04
XIZX = 9.820000e+02
XIXY = 0.000000e+00
XIYZ = 0.000000e+00
XMASS = 6.371595e+02
DYCG = 0.000000e+00
DZCG = 0.000000e+00


def F16_inertia(CG_PCT_MAC):
    DXCG = 0.01*CBAR*(35.0 - CG_PCT_MAC)
    return (numpy.array([XIXX, XIYY, XIZZ, XIZX, XIXY, XIYZ, XMASS, DYCG, DZCG, DXCG]))
